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ABSTRACT 

Experiments were done to learn the critical flowrate of two-phase 
nitrogen through an orifice. Recent analyses of the critical flowrate 
of a two-phase fluid are summarized and compared with experimental 
data. The experimental data appear to fit Ward's analysis better than 
any other analysis. 
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DEFINITION OF SYMBOLS 


Symbol 

Definition 

Units 

A 

Cross-Sectional Area 

m 2 

a, b, c, d, e 

Parameters Defined in Equation (19) 

- 

h 

Enthalpy 

J/G 

K 

Slip Ratio Defined by Equation (13) 

- 

P 

Pressure 

N/cm 2 

Q 

Volumetric Flowrate 

m 3 / s 

R 

Gas Constant 

m 2 /°K kg s 2 

s 

Entropy 

J / G°K 

T 

Temperature 

°K 

V 

Specific Volume 

m 3 /kg 

Vm 

Specific Volume Defined by Equation (Z) 

- 

W 

Weight(Mass) Flowrate 

kg/ s 

X 

Gravimetric Quality 

- 

y 

Void Fraction or Volumetric Quality 

- 

Ah 

Enthalpy of Vaporization 

J/G 

As 

Entropy of Vaporization 

J/G°K 

Y 

Ratio of Specific Heats 

- 

£ 

Weight Flowrate per Unit Cross- 
Sectional Area 

kg/m 2 s 


Fluid Velocity 

m/ s 


Vll 
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Symbol 

g 

l 

m 

o 

P 
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SUBSCRIPTS 

Definition 

Vapor State 
Liquid State 
Maximum 
Stagnation State 
Pressure 
Constant Entropy 
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CRITICAL FLOWRATE OF TWO-PHASE NITROGEN 

By 

Hugh M. Campbell, Jr. and Thomas J. Overcamp* 
George C. Marshall Space Flight Center 


SUMMARY 


Experiments were done to determine the critical flowrate for 
two-phase nitrogen. Recent critical two-phase fluid flow analyses 
are summarized and compared with experimental results obtained in 
this program. 

Most critical two-phase fluid flow analyses utilize an assumption 
of thermodynamic equilibrium, which causes experimental results to 
deviate from the analysis at low qualities. This characteristic devia- 
tion was noted in the data obtained during this program. The experi- 
mental data appear to be more closely approximated by Ward*s 
analysis than any other. 


INTRODUCTION 

Since the end of World War II scientists have devoted much time 
to the investigation of two-phase fluid flow, critical two-phase fluid 
flow, and heat transfer to a two-phase fluid. The major portion of the 
work on two-phase fluid flow was for application to nuclear reactor 
design. Consequently most of the analyses were verified with water 
and steam. With the advent of space boosters, considerable curosity 
has been expressed over the critical two-phase flow for the cryogenic 
fluids. The investigation of two-phase fluid flow phenomena is impor- 
tant in answering questions on the venting of propellant tanks while in 
orbit, chill of transfer and engine propellant lines, and in advancing 
the state-of-the-art in booster design. 

The critical flowrate of a fluid is not a new phenomenon, but is 
well known, both theoretically and experimentally, for a compressible 
fluid. The theory of critical gaseous fluid flow is developed in most 
thermodynamics and compressible fluid-flow textbooks. When a 
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gaseous fluid flows through a restriction, the weight flowrate is 
increased by decreasing the pressure downstream of the restriction, 
because the potential causing the fluid to flow is the pressure gradient. 
As the downstream pressure continues to decrease, the weight flow- 
rate continues to increase until a limiting weight flowrate is attained. 
This limiting weight flowrate is attained when the fluid velocity in the 
throat of the restriction coincides with the speed of sound in the flow- 
ing fluid. Once the critical flowrate is attained, additional decreases 
in downstream pressure will produce no change in the weight flowrate. 
With a constant upstream pressure, the maximum ratio of the pressure 
upstream of the restriction to the pressure at the smallest cross- 
sectional area in the restriction is known as the critical pressure 
ratio. The relationship between weight flowrate and pressure ratio 
is shown in Figure 1. The upstream pressure can be increased when 
the pressure ratio is subcritical, and the only changes in weight flow- 
rate will be due to changes in fluid density or changes in the speed of 
sound in the flowing fluid caused by the increased pressure. 

Two-phase fluid flow exhibits the same critical flow phenomenon 
as gaseous fluid flow. Two-phase critical flow is not analyzed so 
readily as gaseous fluid flow, and consequently, it is not analyzed 
in the textbooks. Critical two-phase fluid flow analysis is greatly 
hindered by thermodynamic metastability. Analyses based on a 
completely homogeneous fluid in thermodynamic equilibrium yield 
results considerably below measured flowrates. Analyses based on 
thermodynamic equilibrium and the various two-phase fluid flow 
patterns yield good results in long tubes where no large amount of 
thermodynamic metastability exists. Analyses assuming no exchange 
of mass from one phase to the other are also used to predict the 
critical two-phase flowrate. 

RECENT ANALYSIS OF CRITICAL TWO-PHASE FLUID FLOW 

Levy [1] predicted the critical two-phase mass flowrate for steam 
and water. Assuming thermodynamic equilibrium, each phase can be 
represented by a unique mean velocity and no frictional or hydrostatic 
head losses. This analysis was based on a momentum exchange 
derived earlier by the same author. 


Using the relation: 



2 



( 2 ) 


The parameter V m is given by: 
v g x 2 vo (1-x) 2 

m y (l-y) 

Using the Bernoulli equation and the equation of continuity with the 
assumptions of no frictional or hydrostatic pressure losses, the 
gravimetric and volumetric quality can be related by: 


X - y u - 2y) + y lu - 2y) 2 + y [(2v g !_vj) (1 - y) 2 + y(l - 2y)j} 1 /2 
[(2 v g /v£ ) (1 - y) 2 + y(l - 2y)] 

Equation (3) is significant because it is an accurate relationship 
between gravimetric and volumetric quality that does not depend on 
the assumption of a homogeneous fluid, and nothing is required to 
relate these two parameters other than the liquid and vapor specific 
volume. 


Equations (2) and (3) were then differentiated to give: 
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These results were then used in the following equation to determine 
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The results of equation (6) were used in equation (1) to determine the 
critical two-phase mass flowrate. This derivation assumes that the 
maximum critical two-phase flowrate is attained when the process 
through the obstruction is isentropic. Comparisons were made for 
an isenthalpic process and the difference between the values obtained 
using these two premises is small. 

Comparison was made between the predicted values and experi- 
mental data from water and steam and the agreement was good. In 
all cases where Levy’s prediction was compared with experimental 
data, the data were slightly below the predicted value or the prediction 
appeared to be a good average of the data. Critical throat pressure 
versus stagnation quality is shown for hydrogen in Figure 2, for 
nitrogen in Figure 4 and for oxygen in Figure 6 using Levy's analysis, 
Critical mass flowrate is shown versus stagnation quality for hydrogen 
in Figure 3, for nitrogen in Figure 5 and for oxygen in Figure 7 using 
Levy 1 s analysis 

Moody [2] predicted the maximum flowrate of a single component 
two-phase mixture assuming annular flow, uniform linear velocities 
for each phase, and thermodynamic equilibrium. He further assumed 
that the liquid and vapor were at the same static pressure at any axial 
location in the nozzle and that the slip ratio (the ratio of the mean vapor 
velocity to the mean liquid velocity) is an independant variable. Using 
the relations: 


i = 


w _ y _ / 1 - y \ h 
A x v g \1 - xj v £ 
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Agr 

V = “A- 


(11) 

1-y - Ai 
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(12) 

k = hA 


(13) 

Equation (7), (8), (9), (10), (1 1), (12) and (13) 
following relationship between x and y: 

can be solved to give the 


y = - K i(4) G*-)l 

By definition, the isentropic stagnation state is defined as, "the state 
a flowing fluid would attain if it underwent a reversibly adiabatic 
deceleration to zero velocity" [3]. Hence: 

S Q = s = s t + x As ( 15 ) 


Equations (7) through (15) were then solved to yield: 
£ = 


(16) 


»0>o - h f > - -££- <s„ - s,)] 


K 2 As 


p + s o- s ^ 


As 


(v g /As) (s o — s t ) 


K v l 
i As 


(s o - 


Equation (16) shows that the mass flowrate through a nozzle is a 
function of stagnation enthalpy, slip ratio, and pressure. Equation (16) 
was differentiated with respect to K and set equal to zero to yield: 

K m =y/ v g /v | (17) 


and was differentiated with respect to pressure, set equal to zero, and 
the following equation resulted: 
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where: 



Curves were then presented showing the local pressure and 
maximum mass flowrates in terms of stagnation enthalpy and stag- 
nation pressure. A comparison was made between the analysis and 
experimental data with a good degree of correlation. Critical throat 
pressure using Moody’s Maximized Analysis is shown versus stagnation 
quality for hydrogen in Figure 8, for nitrogen in Figure 10, and for 
oxygen in Figure 12. Critical mass flowrate using Moody’s Maximized 
Analysis is shown versus stagnation quality for hydrogen in Figure 9, 
for nitrogen in Figure 11, and for oxygen in Figure 13. 



Ward [4] derived two models for predicting the critical two-phase 
mass flowrate. The following equations were used in both analyses. 
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t, = 
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Equations (20) through (23) were combined to yield: 


( 1 - x ) _jc_ 

It ^ 


(24) 


The following assumptions were utilized in the "Frozen Equili- 
brium Model": 

1. No frictional pressure losses 

2. No mass or energy exchange between the two phases 

3. Isentropic expansion 

4. Uniform pressures across any cross section 

5. No shearing stress between the liquid and vapor phase 

6. Negligible velocity at the orifice or nozzle inlet 

7. The vapor phase behaves as a perfect gas 

The liquid and vapor velocities may then be expressed as: 

= \/l0 1 2 3 4 5 6 7 (p o - p) v^ (25) 
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Equations (25) and (26) can then be combined with the continuity 
equation to yield: 



( 27 ) 


(28) 


Equations (27) and (28) were substituted into equation (24), differentiated 
with respect to (p / p^) , and equated to zero to yield: 


d(p/p o ) 



(29) 


-<VPq>1" 5 (¥) Ifefr) <W 1 2 - V)/V - <P,/Po> 1/V 


55 0 


(^)°- ! |<P t /Po» 2 /Y -(P t /Po) (Y+U/Y|, - S 


The ratio of throat pressure to inlet pressure that would satisfy 
equation (29) was then determined for several values of the specific 
heat ratio and inlet pressure. These values of the ratio of throat 
pressure to inlet pressure were then used in equations (24), (27) and 
(28) to compute the critical two-phase mass flowrate. Critical throat 
pressure versus stagnation quality using Ward ! s Frozen Equilibrium 
Model is shown for hydrogen in Figure 14, for nitrogen in Figure 16 
and for oxygen in Figure 18. Critical mass flowrate versus stagna- 
tion quality is shown for hydrogen in Figure 15, for nitrogen in 
Figure 17 and for oxygen in Figure 19. 


The following assumptions were used with the Separate-Phase 
Shifting Equilibrium Model: 


1. No xrictional pressure losses 

2. No mass or energy exchange between the two phases 
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3. Thermodynamic equilibrium 

4. Isentropic expansion 

5. The portion of the liquid which vaporizes or gas that condenses 
continues to move at the same velocity as its parent phase. 

6. The static pressure is uniform across any cross section 

7. Negligible velocity at the nozzle or orifice inlet 

From the general energy equation, the liquid and gas velocity can 
be expressed as: 


+* = 

y 2 x io 4 (h io - h f ) 

(30) 

V 

y 2 x 10 4 (h go - h g ) 

(31) 


Combining equations (30) and (31) with the continuity equation and 
substituting this result into equation (24) yields: 


c = 


(1 - 


x i v c 


y 2 X 10 4 (h io - h s ) / 2 X 10 4 (h go - h g ) 


(32) 


Assuming a constant inlet pressure and quality, then allowing the 
downstream pressure to decrease, equation (32) will attain a maximum 
and thereafter will decrease. This maximum point is the maximum 
flowrate obtained by the "Separate-Phase Shifting Equilibrium Model." 


Critical throat pressure versus stagnation quality using this analy- 
sis is shown for hydrogen on Figure 20, for nitrogen on Figure 22, and 
for oxygen on Figure 24. Critical mass flowrate versus stagnation 
quality is shown on Figure 21 for hydrogen, for nitrogen on Figure 23, 
and for oxygen on Figure 25. 


USE OF CURVES 


Using Figure 2 through 25, the critical pressure and mass 
flowrate may be read directly. The charts using the desired model 
and fluid are first determined. The stagnation quality and pressure 
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are entered, and the critical throat or receiver pressure is determined. 
Any absolute receiver pressure below the critical pressure will not 
alter the mass flowrate. Any absolute receiver pressure above the 
critical pressure will result in sub-critical fluid flow, and the flowrate 
can be determined using various equations for two-phase fluid flow 
through a restriction. 

The stagnation quality and pressure are then entered into the 
charts to determine the critical mass flowrate. The critical mass 
flowrate will be unique when the receiver pressure is below or equal 
to the critical pressure. 


EXPERIMENTAL APPARATUS AND PROCEDURE 


Experimental Apparatus 

The experimental apparatus consisted of a 4. 7 cm (1.87 in.) 
diameter by 86.5 cm (34 in.) long horizontal, cylindrical plenum 
chamber. The two-phase fluid entered the plenum through a 1 . 06 cm 
(0.42 in.) diameter square-edge orifice and exited through a 2.87 cm 
(1.13 in.) diameter square-edge orifice. 

Two-phase fluid flow was generated upstream of the plenum by 
mixing liquid and gaseous nitrogen. Slightly subcooled liquid nitrogen 
was supplied to the region immediately upstream of the plenum from a 
liquid nitrogen Dewar. Liquid flowrate was measured with a turbine- 
type flowmeter. Liquid nitrogen temperature was measured with a 
copper-constanta n thermocouple . Gaseous nitrogen from a high pres- 
sure GN 2 supply system was supplied to the region immediately up- 
stream of the plenum where it was mixed with LN 2 to give the two-phase 
fluid with the desired quality. GN 2 flowrate was measured with a 
turbine-type flowmeter, temperature was measured with a copper- 
constantan thermocouple, and pressure was measured with a pressure 
gauge. The plenum chamber and mixing region were insulated with 
polyurethane foam to prevent heat transfer to the fluid. 

The two-phase fluid was exhausted from the plenum chamber 
through the 2.87 cm diameter orifice into a vacuum chamber. 

A schematic of the experimental apparatus is shown in Figure 26, 
and the apparatus is pictured in Figure 27. 
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Experimental Procedure 

Upon starting an experiment, liquid and gaseous nitrogen were 
mixed in such proportions as to give the desired quality. Fluid was 
then allowed to flow through the mixing region and plenum chamber 
into the vacuum chamber until the apparatus had reached thermal 
equilibrium. After sufficient time was allowed for the apparatus to 
reach thermal equilibrium, the data were recorded. 

To insure critical flow through the upstream orifice, the plenum 
chamber pressure was maintained below 2. 85 N/ cm 2 (4. 0 psia). 
Experimental conditions were selected so that the entire quality range 
was covered. 


Data Reduction 


Volumetric flowrates were converted to mass flowrates using fluid 
properties tabulated in Reference 5 . Knowing the mass flowrates and 
fluid states, and assuming: 

1. Negligible fluid kinetic energy 

2. No heat transfer to the fluid 

the fluid quality can be computed by the following equation. 

x = ^v^vo + ^l^lo _ kjo (33) 

Ah 0 (W V0 + W io ) Ah 0 

Results 

Thirty- six experiments were conducted that covered the entire 
quality range. Experimental results are tabulated in Table I and are 
plotted versus stagnation quality in Figure 28. 

Experimental data are compared to the predictions in Figure 28. 
The experimental data were divided by 0 o 6 to account for the effect of 
the orifice coefficient. 
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CONCLUSIONS 


Upon comparing the analytical predictions to the experimental 
data, it can be concluded that the analytical predictions are valid for 
the cryogens. After adjusting the experimental data to compensate 
for the orifice coefficient,^ can be seen that the experimental data are 
best approximated at high quality by Ward’s Separate Phase Shifting 
Equilibrium Model, and at low qualities by Moody’s Analysis. The 
shape of the experimental data, when plotted versus stagnation quality, 
is best approximated by Ward’s Frozen Equilibrium Model. 

It should be noted that the predictions are valid for a medium to 
high-quality fluid, but fail for a low-quality fluid. 
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TABLE I 


EXPERIMENTAL RESULTS 


Experiment 

No. 

Qi 

(m 3 / s) 

Q g 

(m 3 /s) 

gn 2 

Temp (°K) 

W 

(kg/s) 

Quality 

x o 

P 2 

(N/cm 2 ) 

P 3 

(N/cm 2 ) 

2 

(N/cm 2 ) ' 

1 

. 004558 

. 5863 

272. 22 

. 080 

.477 

17. 0Z4 

2.241 

1.469 

Z 

. 001699 

. 6160 

267. 77 

. 043 

.969 

16. 962 

2. 186 

1.469 

3 

. 004360 

. 60923 

268. 33 

. 078 

.511 

17.265 

2.227 

1.0411 

4 

. 002378 

. 6231 

268. 88 

. 052 

.799 

17. 106 

2. 317 

1.731 

5 

. 005747 

. 5475 

267. 77 

. 094 

. 369 

16. 762 

2.227 

1.496 

6 

. 004756 

. 5772 

267. 77 

. 083 

.456 

17. 106 

2.186 

1. 331 

7 

. 003171 

. 5911 

267. 77 

. 062 

. 639 

16.893 

2.255 

1.641 

8 

. 00198Z 

. 60815 

Z65. 55 

. 047 

.877 

17. 106 

2. 358 

1.800 

9 

. 004869 

. 55Z0 

Z66. 66 

. 083 

.431 

17. 0Z4 

Z. 310 

1. 731 

10 

. 004869 

. 5602 

266. 66 

. 083 

.436 

17. 024 

2. 427 

1.531 

11 

. 002180 

. 5738 

Z66. 66 

. 049 

.801 

17. 148 

2.. 358 

1. 827 

1Z 

. 00ZZ93 

. 5716 

267. 77 

. 050 

.777 

16.824 

2. 158 

1. 331 

13 

• 00Z491 

. 5682 

267. 77 

. 052 

. 733 

16. 962 

2.241 

1. 627 

14 

. 002378 

. 5682 

268. 88 

. 051 

.754 

17. 0Z4 

2. 172 

1.282 

15 

. 005747 

. 4960 

270. 00 

. 094 

.337 

16.955 

2.462 

Z. 000 

16 

. 007Z47 

. 6540 

270. 00 

. 117 

.351 

16. 755 

2.400 

1.731 

17 

. 002689 

. 5682 

271. 11 

.055 

.697 

17. 106 

2. 358 

1.972 

18 

. 005152 

.5419 

271.11 

. 086 

.402 

16.755 

2. 344 

2. 000 

19 

. 008238 

. 1829 

272. 22 

. 114 

. 062 

16.893 

2.310 

1.469 

ZO 

. 004360 

.5900 

271.11 

. 077 

.497 

17. 106 

2.372 

1. 813 

Z1 

. 008238 

.4617 

270. 00 

. 1Z4 

.219 

16. 824 

2. 344 

1.400 

ZZ 

. 013815 

. 5648 

270. 00 

. 09Z 

. 391 
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FIG 2 Critical Throat Pressure Versus Stagnation Quality 
for Hydrogen (Levy's Momentum Exchange Analysis) 
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FIG 4 Critical Throat Pressure Versus Stagnation Quality 
for Nitrogen (Levy’s Momentum Exchange Analysis) 
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FIG 5 Critical Mass Flowrate Versus Stagnation Quality for 
Nitrogen (Levy's Momentum Exchange Analysis) 
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FIG 6 


Critical Throat Pressure Versus Stagnation Quality 
for Oxygen (Levy's Momentum Exchange Analysis) 
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FIG 10 Critical Throat Pressure Versus Stagnation Quality 
for Nitrogen (Moody's Maximized Analysis) 
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FIG 12 


Critical Throat Pressure Versus Stagnation Quality 
for Oxygen (Moody's Maximized Analysis) 
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FIG 14 Critical Throat Pressure Versus Stagnation Quality 
for Hydrogen (Ward's Frozen Equilibrium Model) 
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FIC 16 Critical Throat Pressure Versus Stagnation Quality 
for Nitrogen (Ward's Frozen Equilibrium Model) 
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FIG 17 


Critical Mass Flowrate Versus Stagnation Quality 
for Nitrogen (Ward's Frozen Equilibrium Model) 
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FIG 18 Critical Throat Pressure Versus Stagnation Quality 
for Oxygen (Ward’s Frozen Equilibrium Model) 
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FIG 19 Critical Mass Flowrate Versus Stagnation Quality 
for Oxygen (Ward's Frozen Equilibrium Model) 
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FIG 20 Critical Throat Pressure Versus Stagnation Quality 
for Hydrogen (Ward’s Separate-Phase Shifting 
Equilibrium Model) 
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FIG 21 Critical Mass Flowrate Versus Stagnation Quality 
for Hydrogen (Ward's Separate -Phase Shifting 
Equilibrium Model) 
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FIG 23 Critical Mass Flowrate Versus Stagnation Quality 
for Nitrogen (Ward's Separate -Phase Shifting 
Equilibrium Model) 
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